(N 

O 
O 
(N 

£ 
(N 



o 



O 

o 



Diffusion coefficient for reptation of polymers 
with kinematic disorder. 

Richard D. Willmann 

Institut fur Festkorperforschung, Forschungszcntrum Jiilich, 52425 Jiilich, Germany 
email: r.willmann@fz-juelich.de 

1st February 2008 



Abstract 

We give a lower bound on the diffusion coefficient of a polymer 

chain in an entanglement network with kinematic disorder, which is 

I , obtained from an exact calculation in a modified Rubinstein-Duke 

\& [ lattice gas model with periodic boundary conditions. In the limit of 

infinite chain length we show the diffusive motion of the polymer to be 
c— ») . slowed down by kinematic disorder by the same factor as for a single 

particle in a random barrier model. 
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Introduction: Among the basic problems of polymer science is the 
derivation of large scale properties of entangled polymers from micro- 
scopic properties, such as the molecular weight, which is proportional 
to the chain length L. A scaling argument due to de Gennesa pre- 
dicts for the zero-field diffusion constant D(0) of a reptating polymer 
chain that to leading order -D(O) oc L~ z , z = 2. Computer simula- 
tionsa'H and experimentseru showed an effective scaling exponent z of 
2.28-2.4 for the accessible range of polymer length in contrast to the 
pure reptation prediction. Reptation based theories accounting for 
contour length fluctuations (CLF)Q and contraint release (CR)Q (par- 
tially) explain this behaviour and predict for increasing chain length a 
crossover to z = 2. So far, this region is not experimentally accessible!! 
In terms of the Rubinstein-Duke (RD) modelJ&Ei which incorporates 
CLF, it is possible to compute the proportionality constant after the 
crossoverO'El i.e.: limL->oo D(0)L 2 = W/(2d + 1), where W is an ele- 
mentary hopping rate setting the time scale of defect diffusion and d is 
the dimensionality of the system environment. Finite size corrections 
behave as D(0)L 2 - W/(2d + 1) oc L~P, where 1/2 < p < 1. The 
experimental relevance of the model is shown in.EJ'Ej The RD model 
is an effective model neglecting many effects such as self avoidance of 
the chain or the short time Rouse dynamics. Moreover, the entangle- 
ment network as encountered e.g. in gel electrophoresis is idealized as 
being regular and static. Real entanglement networks have a random 
structure whose effects on the motion of the polymer have to be taken 
into account x3 

• Spatial variations of the mobility of the 'defects' of stored length. 

• Locally fluctuating potential energy due to interactions between 
chain and environment. 

• Entropically favourable regions of low entanglement density. 

• Relaxation of the environment (CR). 

As many effects are at interplay, it is experimentally impossible to 
isolate the influence of a single one. However, theoretical considera- 
tions and computer simulations can be used to investigate each effect 
separately. The reviewE-3 treats entropic effects and the occurance of 
'entropic trapping'. Relaxation of the environment is of minor impor- 
tance in gels but is considered important in polymer melts.0 
The scope of this communication is the investigation of the influence 
of kinematic disorder, i.e., disorder reflecting varying defect mobility 
without affecting the equilibrium configuration of the chain. InEj an 
analysis of Monte Carlo data for a polymer with kinematic disorder, 
i.e., spatially varying mobility of defects, is performed, which shows 
reptation dynamics to prevail. Being based on computer simulations 



and thus short chains, only speculations about the limit of infinite 
chain length are possible. It is argued that in this limit the diffusion 
constant still scales as 

D(0)L 2 = c, (1) 

where the constant of proportionality c might be some average of the 
hopping rates W. Using a modified RD model we can give partial 
confirmation to this conjecture by rigorously proving a lower bound 
jjper on £)(0), which yields (HI). Moreover, we explicitly calculate the 
constant c and thus show, how the disorder changes the coefficient. In 
the limit of infinite chain length D per L 2 = 1/ < 1/W >. 
Definition of the model: In the RD model, the entanglement network 
is represented as a cubic lattice, the lattice constant being equal to the 
mean entanglement length. A string of L + 1 'reptons', i.e., sections 
with a length of the lattice constant, represents the polymer. The 
repton dynamics is as follows: 

a) Each cell occupied by the chain must contain at least one repton 
to ensure connectivity of the chain. 

b) End reptons can move to adjacent cells provided rule a) is not 
violated. 

c) Interior reptons can move to cells occupied by the neighbour- 
ing reptons if allowed by a). This ensures the dynamics to be 
reptation. 

Considering kinematic disorder each boundary between cells has 
assigned to it an individual hopping rate for a repton crossing in any 
direction . We assume there to be a € N possible rates W a , each 
occuring with probability f(W a ) throughout the network (Fig. 1). 
We demand that for the distribution f(W a ) the disorder averages 
(1/W) and (1/W 2 ) are finite. 

The RD model is a model for electrophoresis. The electric field E 
points along a body diagonal of the cubes in the lattice and each repton 
carries a charge. By local detailed balance, this modifies the rates 
for reptons crossing cell boundaries by a factor B depending^ on it 
moving along (+) or against (-) the field, where B = exp(£'/2)JlJ The 
configuration of the chain can by rephrased as a one dimensional lattice 
gas model with L sites by considering the links between reptons with 
respect to E. Links between reptons in the same cell are represented 
as '0' (vacancy), those which are oriented along (against) the field 
and across a cell boundary with rate W a assigned to it as particles 
of type 'a' ('—a'). Thus, the chain conformation is represented by 
L 'pseudospins' y\ to y^ (Fig. 1). Rule c) for the repton dynamics 
enforces the lattice gas dynamics to be as for an exclusion process: 
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Figure 1: Two dimensional representation of a network with a polymer chain 
and mapping to the lattice gas model. Arrows show possible moves. 



In the bulk particles of sort 'dba' hop to the left with rate B W a 
and to the right with rate B^ W a , where each site can be occupied 
by at most one particle. The end dynamics in the lattice gas picture 
needs some care: Assuming y\ (y£) to be non zero, the only possible 
move is the retraction of the end repton to the cell occupied by it's 
neighbour (rule a)). This retraction, being an annihilation event in 
the lattice gas picture, happens with the same rate as the respective 
move in the bulk. Assuming y\ (y£) to be zero, the end repton can, 
according to rule b), move to any of the 2d adjacent cells. For half 
of these the move leads to links being along the field direction, the 
other half against it. The probability of the chosen move leading to 
the crossing of a cell boundary with rate W a being assigned to it is 
/(Wq). Thus the move of the repton, being a creation event in the 
lattice gas picture leads to y\ (jjl) changing from '0' to '±a' with rate 
B^ l f{W a )W a d {B ±1 f{W a )W a d). This choice of boundary dynamics 
is on average correct, but neglects the actual local structure of the 
network.^ Yet to define is the motion of the centre of mass position 
x in terms of the lattice gas model: 

• Particle type 'a' moving to the right (left) decreases (increases) 
re by 1/(L+1), as this is equivalent to a repton moving downward 
(upward). As there are L + l reptons each contributes 1/(L + 1) 



to the centre of mass position. 

• Particle type '—a' moving to the right (left) increases (decreases) 

xby 1/(L + 1). 

In the RD-model, discriminating links between reptons along ('a') 
and against ('—a') the field direction, which is an arbitrarily chosen 
direction in space, allows for following the transport of stored length 
along this direction. Thus the zero field diffusion constant along this 
space direction can be calculated Jj'H3 which immediately yields the 
3-dimensional diffusion constant as diffusion at zero field is isotropic. 
This is in contrast to the original Rubinstein model,Q which allows 
only for the calculation of the curvilinear diffusion constant along the 
contour of the tube within the model and requires additional assump- 
tions to relate it to the 3-dimensional diffusion constant. 
Relation of open and periodic system: Calculations proceed along 
analogous lines as in.BH3 The adaption to the disordered system is 
straightforward, details will be presented in a forthcoming publica- 
tion.!^ Using detailed balance, we calculated the stationary state 
P* n (0) at zero field with respect to H open . It is a product mea- 
sure and the probability of finding a configuration y = (y\, ..,i/l) is 
given by 



iW<» - A Hvi) wi* p( Vi) = { ,$$j£\ „ £ i 



l/(2d + 1) for m = 0. 

(2) 



The shape of the chain only depends on the signs of the y,. According 
to (|2[) the probability for yi = is l/(2d + 1), for y, being positive 
(negative) d/{2d+ 1). These probabilities are as for the original RD 
model, implying that our kind of disorder leaves the equilibrium con- 
formation of the chain unaffected, as inJi^l 

For the same bulk dynamics, but periodic boundary conditions, it 
turns out that P*(0) is at zero field also a stationary state. This en- 
ables comparing the influence of the boundary terms of the stochastic 
generators on the diffusion constant as ir£j and we can prove that 
Dl(0) > -^L+i> where -Dl(O) means -D(O) for a lattice gas with L 
sites. D^+l ls tb- e centre of mass diffusion constant for a lattice gas 
of L + 1 sites and periodic boundary conditions, where the centre 
of mass variable x depends on particle moves in the lattice gas as 
was stated above for the open system. In the following we calculate 
^L+i t° leading order and thus provide a lower bound on -D(O). As 
it is undisputed that disorder slows down diffusion compared to the 
naive approximation c =< W >, this is the physically relevant bound. 
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Figure 2: Zero range particles moving in a random barrier energy landscape. 



Zero range picture: Dealing with a periodic system it is convenient 
to use the following alternative point of view: Instead of characteriz- 
ing the system by y = (yi, -, ])£)■> it is possible to use the sequence 
s = (si, .., sm) of the signs of those in, which are non zero and have 
rates w = (w±, ..,wm), Wj £ {W±, .., W a }, and n = (m, ..,um), where 
rij equals the number of yi = between Sj and Sj+i- Here, sm+i = s\. 
rij can be understood as counting particles of zero interaction range 
at site j. To make this alternative picture consistent with the lat- 
tice gas picture, we require X^i=i n j = L — M. Also the dynamics 
of the rij is thus fixed: The configuration (..,rij,rij + i, ..) changes to 
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l,rij + i + l, ..) with rate B s i +1 Wj+i and to (..,nj + l,rij+x — l, 



with rate B~ s i +1 Wj+i. This means that the random hopping rates as 
well as the Sj are not assigned to individual particles, but to bonds be- 
tween sites in the zero range (zr) picture.O At E = the zr-particles 
move as in a random barrier energy landscape (Fig. 2). A zr-particle 
hopping to the right across a bond with Sj > (sj < 0) increases (de- 
creases) the centre of mass position by 1/L. Conversely, a zr-particle 
hopping to the left across a bond with Sj > (sj < 0) decreases 
(increases) the centre of mass position by 1/L. Regarding a periodic 
system, the phase space is non ergodic, as in the zr picture neither the 
number, nor the order of the Sj and Wj can be changed. Therefore, the 
phase space constists of connected subspaces, called 'channels'!?! being 
characterized by s and w. For obtaining the expectation value of the 
centre of mass drift velocity for the full phase space v, at first the 
expectation value of the centre of mass drift velocity for each channel 
< v > has to be calculated. Subsequently averaging over the expecta- 
tion values for each channel yields v. To compute D v ^ r we will employ 
the Einstein relation D v ^ r = 1 / 'L(dv / 'dE) e=o- 

Calculations for individual channels: Expressing the dynamics for the 
zr-particles as stated above by a stochastic generator H z i- , we com- 
puted the stationary state Pzr" for arbitrary E. The use of a product 
measure ansatz leads to a recursion relation yielding the following 



steady state probability for a configuration n = (m,..,n M )E 

M i M 

^ w (n) = IK7£]>?- 0) 

j=l n j=l 

The primed sum means summing with the constraint Y2j=i n j = L—M 
and 

M j-1 

Zj = ^2 l/(exp(s j+i E/2)w j+i ) Y[ exp(-s j+k E). (4) 

t=l fc=i 

The drift velocity is in the lattice gas picture given by the difference of 
currents of particles with yi < and yi > 0: < v > = < j~ —j + >• This 
translates into the zr picture as the current of zr-particles across bonds 
with Sj > minus the one across bonds with Sj < 0. Due to using the 
Einstein relation, only first order terms in an expansion of < v(s, w) > 
into E contribute to D p ^ r . The quantity Zi as occuring in (||) can then 
be calculated as Zi = Ylj=i ^/ w i+j = z - This facilitates evaluating 
< v(s, w) >, which using the quantum Hamiltonian formalismE3 and 
a treatment as inH yields 

< 1,(8, w) >= ES 2 ~I^L\ + o(E 2 ) (5) 

with S = J2j=i s j- 

Averaging over the channels: To obtain v, an average of < v > over the 
channels has to be performed, where each channel has to be weighted 
such that in zr and lattice gas picture corresponding configurations 
have equal weight in the stationary state, ^(s, w) is the weight factor 
of the channels as in,BE3 which is modified by the disorder to 
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Thus v is to first order in E given by 



Yl Yl Yl *(S,W)<U(S,W)> 

M s=(si,...,s M )w=(wi,...,w M ) /j\ 
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where < 1/z > is the average of z with respect to the distribution 
f(W a ) and v^ is the result for the RD model without disorder. For 
M — > oo the restrictions on f(W a ) allow invoking the central limit 



theorem, which yields: < 1/z >= 1/(M < 1/W >), leading, when 
employing the Einstein relation, to the following result in the limit of 
infinite chain length: 

^^ = (2d + l)<l/W>- (8) 

Comparing to the result ina for the ordered case RD model, D p " L 2 = 
l/(2d + 1) reveals the remarkable result that in this limit the cen- 
ter of mass diffusion in the RD model on a ring is slowed down by 
kinematic disorder in the same manner as the single particle diffusion 
constant in a random barrier model. We remark that for the ordered 
RD model, ^jr,+i anc ^ Dl{0) are equal to leading order as shown intL2l 
by a variational statement for -D(O). This variational techniquejs also 
applicable to the RD model modified by kinematic disorder^ and 
shows, that the naive approximation c =< W > provides an upper 
bound for -D(O). We performed Monte Carlo simulations of the model 
with open boundary conditions with various distributions f(W) and 
chain length up to L = 40. Apart from quickly decaying finite size ef- 
fects, which depend on f(W), the results clearly indicate that also for 
the RD model with kinematic disorder open and periodic system have 
to leading order the same zero field diffusion constant and therefore 
indeed c = 1/ < 1/W > yields the correct asymptotic behaviour. 
Gunter M. Schiitz is gratefully acknowledged for posing the problem 
and many helpful discussions. 
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